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O and early B stars are at the apex of galactic ecology, but in the Milky Way, 
only a minority of them may yet have been identified. We present the results of a 
pilot study to select and parametrise OB star candidates in the Southern Galactic 
plane, down to a limiting magnitude of ^ = 20. A 2 square-degree field capturing the 
Carina Arm around the young massive star cluster, Westerlund 2, is examined. The 
confirmed OB stars in this cluster are used to validate our identification method, based 
on selection from the {u — g^g — r) diagram for the region. Our Markov Chain Monte 
Carlo fitting method combines VPHAS-F ix, g, r, i with published J, i7, K photometry 
in order to derive posterior probability distributions of the stellar parameters log(Teff) 
and distance modulus, together with the reddening parameters Aq and Ry- The stellar 
parameters are sufficient to conhrm OB status while the reddening parameters are 
determined to a precision of cr(Ao) ^ 0.09 and (t{Rv) ^ 0.08. There are 489 objects 
that ht well as new OB candidates, earlier than ^B2. This total includes 74 probable 
massive O stars, 5 likely blue supergiants and 32 reddened subdwarfs. This increases 
the number of previously known and candidate OB stars in the region by nearly a 
factor of 10. Most of the new objects are likely to be at distances between 3 and 6 
kpc. We have conhrmed the results of previous studies that, at these longer distances, 
these sight lines require non-standard reddening laws with 3.5 < Ry < 4. 

Key words: stars: early-type, (Galaxy:) open clusters and associations: individual: 
Westerlund 2, (ISM:) dust, extinction. Galaxy: structure, surveys 


1 INTRODUCTION 


Stars of spectral type O and early B, more massive than 
~ 8 M 0 , are massive enough to form collapsing cores at the 
end of their nuclear-burning lifetimes (see e.g. Langer|2012 


E-mail: m.smithlO@herts.ac.uk 


Smartt 2009). It is widely recognised that these stars - hence¬ 


forward OB stars - are an important source of kinetic energy, 
driving turbulence and mixing of the interstellar medium, 
powered by a range of phenomena (stellar winds, wind-blown 
bubbles, expanding HII regions and supernova explosions). 
They are the main source of ultra-violet radiation in galaxies 
and, being short-lived (< 40 Myr), they are excellent tracers 
of recent star formation. 
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In the Galaxy, clusters containing OB stars and OB as¬ 
sociations have played an important role in tracing spiral 
arm structure (e.g. Russeil|[20Q3 Vallee|[2008 ). The typical 
scale height estimated for OB stars, forming in the Galac¬ 
tic disk, is a few 10s of pc (e.g. Reed|[2000 Garmany et al. 


Drew et al.|2014 ). VPHAS+ is a deep, uniform, photometric 


1982 ), in keeping with estimates of the scale height for giant 
molecular clouds, their birth sites (e.g. Stark Lee||2005 ). 
OB stars are usually regarded as forming in clustered en¬ 
vironments ( Zinnecker Yorke||2007 ) and are less common 
in the field. However, examples of isolated field O stars are 
known and the question has arisen as to whether these high- 
mass stars have formed in situ, perhaps as the result of 
stochastic sampling of the initial mass function (IMF) as 
outlined by Parker Goodwin] (2007), or have been ejected 
from clusters as runaways (see e.g. Portegies Zwart et al. 
2010 Bestenlehner et ar]2011| . In the Milky Way ~ 96% 
of known 0-type stars have been identified as members of 
young open clusters, OB associations or as otherwise kine¬ 


matically linked to clustered environments (de Wit et al. 


2005). This leaves up to ~ 4% of Galactic 0-type stars pos¬ 
sibly forming in isolation. Deep comprehensive searches for 
OB stars away from clusters have not been undertaken hith¬ 
erto. 

As luminous objects detected to great distances across 
the Galactic disk and through substantial obscuration, OB 
stars have long been recognised as a highly-suitable means 
for characterising the spatial variation of interstellar extinc¬ 
tion, in terms of both dust column and extinction law (e.g. 
Gardelli et al. 1989 Fitzpatrick &: Massa|2QQ7| ). This is aided 
by their relatively simple optical near-infrared (OnIR) spec¬ 
tral energy distributions (SEDs). It follows from this that 
the more densely we can map the positions and extinctions 
towards these luminous probes, the more high-quality em¬ 
pirical constraints we can set on the 3-D distribution of dust 
and dust properties across the Galactic Plane. 

Both of the above areas of enquiry will be well served 
by a deeper, more comprehensive mapping of the OB stars 
in the Milky Way. Past cataloguing efforts have been limited 
to brighter, nearer objects (e.g Garmany et al.p982| Reed 


2003 Mafz-Apellaniz et al.|20Q4 ). Indeed the most compre¬ 

hensive collection so far, ‘The Gatalog of Galactic OB Stars’ 
( Reed|20Q3 ) contains ~ 16000 known or suspected OB stars 
taken from across the literature: around 95% of the entries 
are brighter than 13^^ magnitude in the visual bands. Now 
is the right time to push the magnitude limit much fainter, 
to ~ 20^^ magnitude, given the likely delivery of astrome¬ 
try to this depth by the Gaia mission from ~ 2017 onwards 
(both parallaxes and proper motions, for details on expected 


performance see de Bruijne 2012). Efficient, purely photo¬ 
metric selection of OB stars in the field as well as in clusters 
continues to be best undertaken at blue optical wavelengths. 


where colour selection via the Q method (initiated by John- 
son Morgan|p~953a ) is proven to separate O and early B 
stars from later type stars. 

The practical motivation of this paper is to establish a 
method of photometric selection and analysis that can form 
the basis for a new homogeneous census of Galactic OB stars 
as faint as ^ ~ 20. Based on a restrained extrapolation of 
the first results presented here, we can surmise that a new 
census will more than double the numbers known. A suitable 
source for the new census will be the VST Photometric Ha 
Survey of the Southern Galactic Plane and Bulge (VPHAS+ 


survey of the entire southern Galactic Plane and Bulge in 
broad-band g, r, i and narrow-band Ha filters on ESO’s 
VLT Survey Telescope (VST). The survey footprint includes 
the entire southern Galactic Plane within the Galactic lati¬ 
tude range of |5| < 5°. The VST’s OmegaGam imager pro¬ 
vides a full square degree field of view with very good spatial 
resolution (0.2” pixels sample a median seeing of 0.8 - 1.0 
arcsec in the ujgjr bands). 

Here, we present a first study that uses broadband 
VPHAS+ data to select and parametrize OB stars in a ~2 
square-degree area, roughly centred on ^ = 284°, b = —0.7°, 
in the part of the Plane containing the young massive clus¬ 
ter, Westerlund 2 (Wd 2), the larger associated HII region 
RGW 49, and the diffuse nebula NGG3199 (see Eigure^. 
Previous optical and near-infrared studies on the stellar con¬ 
tent of Westerlund 2 have focused on the immediate environ- 


(Moffat et al. 1991 

Ascenso et al. 2007 

Vargas Alvarez et al. 

2013), while the x- 

-ray study by Tsujimoto et al. (2007 

) fo- 


cused on an area ~ 17 arcmin across. Most recently. 


et al. (2015) have revisited optical photometry of this clus¬ 


ter over a 17.9’ x 9.3’ footprint. 

By tracing Sgm warm-dust emission [Rahman V Mur^ 
ray (2010) have identified this same region as part of a large 


star-forming complex (G283). On the sky, Wd 2 falls close 
to the Garina Arm tangent direction (e.g. Russeil|2003 ): the 


Dame 


(2007) show persuasively that 


GO data presented by 
Wd 2 and its environs fall just inside the sky position of the 
tangent point, but fu rther away. This cluster is estimated 


to be 1 - 3Myr old (Vargas Alvarez et al. 


2013 


Ascenso 


et al.|2007 ). It contains a large number of spectroscopically- 


confirmed OB stars, albeit behind a dust column giving rise 
to over 6 magnitudes of visual extinction ([Moffat et ajjl99f 


Rauw et al.|2007 Garraro et al.|2012 Vargas Alvarez et al. 

2013| ). Estimates of the distance to Wd 2 in the literature 


have varied enormously, ranging from 2.8 kpc ( Ascenso et al 
2007) up to ~8 kpc (e.g. Rauw et al.||2011 ). However, it is 


not our aim to enter into this debate. More important is the 
likelihood that much of the scientific gain from VPHAS+ 
discoveries of OB stars will be in the domain of visual ex¬ 
tinctions of up to 8-10 magnitudes, and distance scales of 
2-10 kpc (according to Galactic longitude). In this regard, 
the field around Wd 2 is highly typi cal of the task ahead. 


A recent study on Wd 2 by Vargas Alvarez et al. 


(2013) uses data from the Hubble Space Telescope (HST) 
that offers much better spatial resolution than is achiev¬ 
able from the ground. This is the only dataset that offers 
better angular resolution than the new VPHAS+ data anal¬ 
ysed here. These authors’ values of Rv and Ay were de¬ 
rived by fitting, to 32 individual OB stars in or near Wd 
2 , reddened model optical/near-infrared SEDs appropriate 
for the selected stars’ spectroscopically-confirmed spectral 
types. The best fits were computed by seeking the global 
chi-squared minimum among all plausible values of Rv and 
Ay - resulting in a mean outcome of Ry = 3.77 ± 0.09 and 
Ay = 6.51 ± 0.38 mag combining results from different red¬ 
dening law prescriptions. We use a comparison of our OnIR 
SED fit results for this same set of OB stars to bench-mark 
our method. 

This paper is organised as follows. In Section 2 more 
details on the data used for this study are given. Section 3 
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Figure 1. RGB image of the ~2 square degree region (Ha, g, i). This region falls within the star forming complex G283 identified by 
[Rahman &: Murray] ( |2Q10| > - an elliptical region slightly larger than the sky area shown. Westerlund 2 is embedded in the HII region 
RGW49, while the diffuse nebulae NGG3199 is located to the right (West) as marked. The dashed line traces the Galactic equator. 


is a presentation of our method, beginning with the updated 
version of the Q method of OB star selection that we use, 
and ending with a description of the Markov Chain Monte 
Carlo (MCMC) sampling of the posterior distributions of 
the OnIR SED model fit parameters. The stage is then set 
to compare our resu lts for Wd 2 stars with those of ] Varga’s] 
Alvarez et al. (2013), in Section 4. The results of the fits to 


the final list of 527 new OB candidates drawn from across 
the full 2 square degrees are presented in Section 5. This is 
followed by a discussion of the results in Section 6, in which 
we consider the extinction trends revealed in this region, 
and draw attention to the newly discovered O stars outside 
the confines of Wd 2. The outlook and our conclusions are 
summarised in Section 7. 


2 THE DATA 

We make use of the photometry from two VPHAS+ fields, 
numbered 1678 and 1679, that are respectively centred on 
RA 10 18 10.91, Dec -58 03 52.3 (J2000) and on RA 10 
25 27.27, Dec -58 03 52.3 (J2000). These were observed in 
succession in the g and r filters on the night of 22^^^ 
January 2012. The red filter data in Ho, r and i were ob¬ 
tained on 29^^ April 2012. The seeing, as measured from 
the data point spread function, was variable on the earlier 
night ranging from 0.62 at best in g up to 1.24 at worst 
in r. When the exposures in the red filters were obtained 3 
months later, conditions were more stable, with the typical 


seeing ranging from 0.8 to 1.0 arcsec. Viewed in comparison 
to all the VPHAS+ data collected so far, these observations 
rank as 2’^'^-quartile quality in u and g (i.e. relatively high 
quality), and S’^^-quartile in r, z and Ho. The 5cr magnitude 
limits on the single exposures are u\ 21.0, g\ 22.4, r: 21.5, i\ 
20.6, and Ho: 20.4. All magnitudes are in the Vega system. 
Full details on the survey strategy, the offsets, the exposure 
times, photometric quality and the data-processing pipeline 
used are given by Drew et al. (2014). 


Our analysis begins with band-merged catalogues cre¬ 
ated from the single-band catalogues emerging from the 
CASH pipeline. In order to correct for the uncertainty in 
the initial calibration of VPHAS+, a comparison has been 
made with empirical r and i observations from the APASS 
survey and with synthetic tracks in the {u — g^g — r) plane. 
The median difference between r and i in the two sur¬ 
veys was applied to the VPHAS+ data. The u band was 
then calibrated by applying an offset to the u — g scale such 
that the number density of stars between the synthetic GOV 
reddening track and the unreddened main sequence is max¬ 
imised. This ensures that the top and bottom edge of the 
main stellar locus are aligned with the synthetic tracks as 
shown in Figure This resulted in offsets relative to the 
pipeline reduction of u\ -0.35, g: 0.05, r: 0.01 and i: 0.05 for 
field 1678 and u: -0.34, g: 0.06, r: 0.01 and i: 0.01 for field 
1679. With an improved calibration in place, we select stars 
in the magnitude range 13 < ^ < 20 and require random 
photometric errors to be less than 0.1. Mean magnitudes 
were taken when repeat photometry was available from the 
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Figure 2. Calibration of VPHAS+ data with respect to synthetic 
reddening tracks from |Drew et al.| ( |2014| ). Both the main sequence 
and GOV reddening vector line up with the main stellar locus. 


offset fields. Objects were removed if the photometry in the 
offset field differed by > 0.2 mags. This removes unreliable 
photometry due to objects that fall on a CCD edge. 


3 SELECTION AND FITTING METHOD 

3.1 Photometric selection and cross matching 


We select OB stars using a method that has its origins in 
the Q Method of Johnson & Morgan (1953b). On the {u — 
g,g — r) diagram reddened OB stars of spectral type earlier 
than B3 are located above and away from the main stellar 
locus. We initially select our candidate objects above the 
reddening vector associated with a B3V. In principle no star 
can be bluer than the Rayleigh-Jeans (RJ) limit which sets 
an upper bound on the likely location of OB stars in the 
diagram. The blue objects that lie above the RJ reddening 
vector were nevertheless included in the selection and their 
origins are discussed in section 

Figure shows the selection of OB candidates (blue 
cross es) across the two fie l ds as well as the known OB stars 


from Vargas Alvarez et al. (2013) that were successfully cross 
matched with VPHAS+ (shown as red triangles). Over¬ 
plotted are the reddening tracks of a B3V, a BIV and that 
of a pure RJ spectrum all taken from [Drew et al.| ( |2Q14 ). 
The tracks we use take into account the measured red leak 
associated with the i^-band filter. 


Previous r esults from Carraro et al. (2012) and Vargas 


Alvarez et al. (2013) suggest an Rv = 3.8 reddening law is 
required towards Wd 2. The BIV and RJ reddening vectors 
have been drawn using this law. To avoid a bias towards this 
non-standard reddening law we have used the B3V Rv = 

3.1 reddening vector as our lower selection limit and have 
dropped its position by 0.1 mags in —^ in order to capture 


Figure 3. Selection of OB stars in and around Wd2. The lower 
reddening curve is that of a B3V, dropped by 0.1 in u-g in order to 
capture all early type B stars, and is characterised by an Ry = 3.1 
law. The other reddening curves are that of a BIV and an ideal 
Rayleigh-Jeans spectrum and are characterised by an Ry = 3.8 
law. Selecte d OB candidates are blue c rosses while the known 

| |2013t 


objects from Vargas Alvarez et al. 


any early B stars that may have been missed. The lower the 
value of Ry, the steeper the reddening vector will be. 

Each object was then cross matched to within l" of 
the best available near infra-red detection in order to access 
J, iJ, K photometry. The mean angular cross-match distance 
was 0.09”. As the stellar density in the central ~ 4' of Wd 
2 is very high, the Ascenso et al. (2007) NIR catalogue was 
the preferred partner on account of its superior angular res¬ 
olution. Everywhe re else 2MASS was used. Thi s follows the 
approach taken by Vargas Alvarez et al. (2013). 


3.2 SED fitting 

We calculate the probability distribution of a range of model 
parameters corresponding to a set of empirical measure¬ 
ments, in a Bayesian scheme. This approach is chosen over 
a straight forward minimisation scheme so that we may 
recover the full posterior probability distribution. This can 
reveal covariance between different parameters. 

Given a set of empirical data, d = {di,..,di}, and a 
model, parametrised by a set of parameters, 0 = {^i, ..,^z}, 
the posterior probability of the parameters can be calculated 
using Bayes’ Theorem: 


P{0\d) 


P{d\0)-P (0) 

W) 


( 1 ) 


In this expression, P (d \ 0), the likelihood is the prob¬ 
ability of the data being measured given a set of model pa¬ 
rameters. The posterior and the likelihood are related by the 
prior, P{0), which encodes any known constraints on the 
model parameters, including known physical bounds. Here 
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P(d) can be treated as a normalising constant and ignored. 
Hence the posterior probability distribution can be found 
by the relation: 


P(0 \d) o^P{d\0)‘P{0) (2) 

In this work the empirical data are derived from the 
observed SED of each star and they consist of optical and 
near infrared apparent magnitudes: 

SEDobs = {u, g, r, z, J, H, Ks}, (3) 

and their uncertainties: 


Cfi^SEDohs) — ■> (^r-) J1 •) • 


(4) 


Along with the random flux errors supplied by the sur¬ 
veys, we have included a systematic uncertainty to account 
for the independent absolute calibration errors in each band. 
The values adopted for the latter are 0.04 in the u band, 0.03 
in g,r and z, 0.03 in the J band and 0.02 in H and Kg (see 
Drew et al.||2014 Skrutskie et al.||2006 ). 


The model parameters that we are interested in esti¬ 
mating are: 


0 = {log(Teff), Ao, Rv, g} (5) 

Where log(Teff) is the effective temperature, Aq is the 
monochromatic extinction at 4595A, Rv is the ratio of total 
to selective extinction and g is the distance modulus. 


Table 1. Sample values of the intrinsic SEDs with approximate 
spectral type equivalents. Magnitudes are in the Vega system. 


ST 

log(Teff) 

u 

9 

r 

i 

J 

H 

Ks 

03V 

4.65 

-7.32 

-5.78 

-5.48 

-5.33 

-4.88 

-4.73 

-4.63 

09V 

4.50 

-5.28 

-3.86 

-3.60 

-3.45 

-3.03 

-2.90 

-2.80 

BIV 

4.40 

-3.97 

-2.70 

-2.47 

-2.34 

-1.98 

-1.85 

-1.77 

B3V 

4.27 

-2.31 

-1.33 

-1.16 

-1.07 

-0.80 

-0.70 

-0.65 


effective temperature (Martins et al. ]200^ , as the Rayleigh- 
Jeans limit is approached. This means that log(Teff) is only 
weakly constrained, albeit well enough to reach our goal of 
confirming OB status. As we have no handle on luminosity 
class, the distance modulus takes the role of a normalisation 
factor and will also be weakly constrained. In contrast Aq 
and Rv are very informative and well constrained. 

We can now use the forward model to construct a like¬ 
lihood model P {SEDobs \ 0) that computes the probability 
of SEDobs given the set of physical parameters 6. Assuming 
that the uncertainties on the measurements are normally 
distributed and uncorrelated, this can be described by a 
multi-variate Gaussian: 


P {SEDobs I 0) oc 


exp 


{SEDobs - SED^odf {SEDobs - SED„ 


( 6 ) 


Where E is the covariance matrix containing the vari¬ 
ance {SEDobs) in the leading diagonal. In this case Equa¬ 
tion reduces to the familiar sum for 


3.2.1 Likelihood funetion 

Defining the likelihood function requires us to define a for¬ 
ward model SEDmod{0), which predicts the apparent SED 
of OB stars based on the model parameters 0. The intrinsic 
SEDs used in the model are taken from the Pa dova isochrone 
database (CMD v2.2 0 [Bressan et al.|[^012 


Bertelli et al. 


1994) and are supplied in the Vega system. The optical/NIR 


colours of OB stars do not vary significantly with luminosity 
class ( Martins et al.|[2QQ5 ). Therefore log(^) was fixed and 
only main-sequence models were used {\og{g) ~ 4.0). Solar 
metallicity Z = 0.019 has been adopted throughout, in view 
of the fact that the sight lines we explore do not sample 
a wide range of Galactic rad i i. Thi s is the same value as 
used by Vargas Alvarez et al. (2013). Eixing these parame¬ 


ters provides a simple grid of absolute magnitude, Ma, as a 
function of log(Teff) in each of the seven bands. 

To obtain a continuous grid, each Mx — log(Teff) re¬ 
lationship was fit with a 2'^^ order polynomial. It can be 
noted that a linear fit was also trialled but failed to charac¬ 
terize the distributions especially for the low-end values of 
log(Teff). Table provides sample SEDs. 

The SEDs are then reddened using a |Eitzpatrick &| 
|m assa| ( [2QQ7| reddening law, parametrised by Aq and Rv, 
and then shifted according to a distance modulus. The ap¬ 
parent OnIR SEDs of O and early B stars are largely con¬ 
trolled by these quantities. This is because the OnIR intrin¬ 
sic colours of OB stars change very slowly as a function of 


^ http://stev.oapd.inaf.it/cgi-bin/cmd 


/ n 

P {SEDobs I 0) ocexp 


{m{obs)i — m{mod)i) 


(7) 

Where m{obs)i and m{mod)i are the observed and 
model magnitudes in each band z. 


3.2.2 Priors 

We adopt a uniform prior on each of the model parameters: 


P{0) 


i 4.2 ^ log(Teff) ^ 4.7 
0 ^ Ao ^ 15 

2.1 ^ Rv ^ 5.1 (8) 

0 ^ g ^ 20 

^ 0 else 


The upper bound on log(Teff) is governed by the avail¬ 
able models and the lower bound is slightly less than the the 
typical temperature of a B3V star ( Zorec BriotlpT^ ) in 
accordance with our selection in the {u — g,g — r) diagram. 
The constraints on Rv are the upper and lower limits mea¬ 
sured in the Galaxy ( Eitzpatrick Massa|2QQ7| . The upper 
limit on Ao is much larger than maximum extinction allowed 
for detection of OB stars in VPHAS+ down to g = 20, as¬ 
suming a typical rise in visual extinction of 1 magnitude per 
kpc. This makes the prior on Ao essentially unbound. The 
upper limit on the distance modulus g of 20 is well beyond 
the realms of the galaxy and so is also essentially unbound. 
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Placing a large but finite limit on Aq and fi enables the 
MCMC algorithm to converge more quickly. 


3.2.3 Sampling the posterior distribution using MCMC 


Characterising the posterior distribution by computing the 
probability at all values in the parameter space is computa¬ 
tionally expensive. Instead one can sample the distribution 
using an MCMC algorithm. 

In this study we use the Python package emeee devel¬ 


oped by Foreman-Mackey et al. (2013). In brief, the software 


takes a set of parameters and supplies them to a group of 
n walkers. The walkers then use a pseudo-random walk to 
sample the parameter space. At each sample the probabil¬ 
ity is calculated. By communicating their relative probabil¬ 
ities to one another the walkers are able to quickly find and 
sample the region of high probability without wasting com¬ 
putational time on the parameter combinations of very low 
probability. The software then returns what are known as 
ehains which contain the values of the parameters at every 
step in the walk. The frequency at which each region in the 
parameter space is visited is proportional to its probabil¬ 
ity. The finer details can be found in|Foreman-Mackey et al.| 
(20131. 


4 VALIDATION OF METHOD 

First it is appropriate to verify that our selection method 
recovers known objects. Second we verify that the fitting 
algorithm delivers the expected results. To achieve this, we 
ha ve chosen to compare with t he results of the recent study 
by Vargas Alvarez et al. (2013). This is an informative com¬ 


parison to make both because this study benefited from 
the superior angular res olution of HST and because |Var-| 
gas Alvarez et al. (2013) used a combination of optical and 
NIR photometry to derive stellar reddenings as we do here. 


4.1 Photometric selection 


Vargas Alvarez et al. (2013) derived the extinction proper¬ 


ties of 29 known OB stars in the central region of Wd 2, 
of which, 24 were successfully cross matched with VPHAS+ 
to within 1”. Using the nomenclature from Vargas Alvarez 


et al. (2013), the five missing objects are #597, #826, #843, 
#903 and #906. They appear in some of the most crowded 
regions of the cluster: the angular resolution of VPHAS+ 
compared to that of HST is insufficient to separate them 
from brighter neighbours. Figure ^ shows the positions of 
the 24 cross-matched o bjects and the positions of those that 
are missing (relative to Vargas Alvarez et al.|2013| over plot¬ 
ted on the g-band image. 

Figure is the highly magnified section of Figure that 
contains the objects with known spectral type. The red and 
blue shaded regions are where we expect to find late-type 
(09 - 06) and early-type (06 - RJ) O stars respectively. 
We find that the majority of the objects are correctly sepa¬ 
rated into their respective early or late spectral-type zones 
defined by the Rv = 3.8 reddening tracks. This gives an 
early indication that an Rv ~ 3.8 reddening law is required 
for this sight-line and that the calibration of the data is in 
good agreement with the synthetic photometry. 



Figure 4. Inverse VPHAS+ g band image of the central region 
of Wd 2 showing the obj ects with known spectral type from | VaiT] 


3ie( 

i- 


gas Alvarez et al. ( 2013|. The red triangles are the positions of 


the objects detected in VPHAS+ and the blue squares are the 
positions of those that are not detected due to crowding. 


Object #771 falls well above the ‘RJ limit’. As a con¬ 
firmed 08V star, its position in the (u — g^g — r) diagram 
is clearly anomalous. Close inspection of the image suggests 
that the photometry of this star is affected by a bright neigh¬ 
bour. 


4.2 SED fitting 

Ultimately 21 of the 24 known objects were suitable for SED 
fitting. These objects are tabulated in Table Two of the 
objects left out are #896 and #771 for which there is no de¬ 
tection in one or more of the the optical bands due to blend¬ 
ing. The third is object #1004 for which the near-infrared 
photometry is incomplete. 

For each of the 21 objects for which we have computed 
SED fits, the posterior distribution was sampled with 100 
walkers over 10000 iterations with a 1000 iteration burn in. 
The typical autocorrelation time for each walk (or number 
of steps per independent sample) was found to be well be¬ 
low 100, which indicates that the posteriors are thoroughly 
sampled. We can determine the probability distributions for 
each parameter by marginalising over all other parameters. 
We visualize this by constructing 1-D histograms of the val¬ 
ues of each parameter visited in the random walk. We can 
also check for covariance or degeneracy between parameters 
by constructing marginalised 2-D histograms for each pair 
of parameters. Eigurej^ shows an example of these diagrams 
for an 04V and a BIV star in the sample (#913 and #549). 

The obvious difference between the two cases is appar¬ 
ent in the 1-D marginalisation of parameters. We see that 
the hotter the object the more skewed the probability distri¬ 
butions in log(Teff) and /a become. This can be attributed to 
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Figure 5. Testing the selection process of OB stars associated 
with Wd2. Objects with known spectral type tend to fall into the 
correct synthetic spectral type range with an Ry = 3.8 reddening 
law. 


the fact that the hotter SEDs are approaching the RJ tail. 
This makes it more difficult to differentiate the temperature 
of the hottest stars and consequently the luminosity and dis¬ 
tance. This makes the drop off in probability at the hot end 
more shallow. This intrinsic feature also means that the un¬ 
certainties on log(Teff) and /j. increase with temperature but 
has the positive effect of decreasing the uncertainties on Aq 
and Rv- For the later type stars log(Teff) is better defined 
but still uncertain. 

The value adopted for each parameter is the median 
of the marginalised posterior distribution with upper and 
lower uncertainties defined by the 16^^ and 84^^ percentiles. 
We find that we are able to determine the values of Aq and 
Rv with relatively high precision (better than ±0.09 mag 
and ±0.08 respectively in all cases). These uncertaintie s are 
similar to those found by Vargas Alvarez et al. (2013). We 
note that Rv and Aq are well defined and show negligible 
covariance relative to each other and only modest covariance 
with respect to log(Teff) and //. 

However, as expected, our determination of tempera¬ 
ture and distance are not so informative. For object #913, 
log(Teff) = 4.611 q;q 7 and = 13.23lQ;8g. This corresponds 
to values of Tgff = 40.7lg ?kK, or a spectral type range from 
08V to 02V. The results for fi translate to d = d.dij'gkpc. 
This already significant distance uncertainty is nevertheless 
an underestimate given that neither the luminosity class or 
metallicity uncertainties have been formally incorporated. 
In addition we are treating all stars as if single which biases 
the inferred distance moduli to lower values by up to 0.75 
mag. Because of the relative lack of constraint on log (Tgff) 
from the intrinsic colours of OB stars, the error in log(Teff) 
is driven mainly by the error in fi. In comparison the direct 
effect of binarity on log(Teff), through colour-changes, will 
be small. It is plainly apparent in Figure that log (Tgff) 








3 , 


log{Teff)[K] = 


4.61 


+ 0.06 

- 0.07 



log{Teff)[K\ A, 


Rv 





logiTeff][K] ^0 Rv n 


Figure 6. PDFs of the fitting parameters as a result of the 
MCMC simulation for stars #913 an 04V ftop) and #549 and 


BIV (bottom) using the numbering system from Vargas Alvarez 

let al.| ( |207^ . 


and fi are strongly and positively covariant. The role of the 
distance modulus is essentially that of a normalisation pa¬ 
rameter. 

Figure shows the results for the 04V star from Figure 
[^translated into the original SED data space. The top panel 
shows the observed SED over-plotted by 30 randomly sam¬ 
pled model SEDs that are drawn from the posterior distribu¬ 
tions shown in Figure The lower panel shows the residuals 
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Wavelength 


Figure 7. Visualisatio n of the posterior distributi ons of objects 
^913 and #549 (from Vargas Alvarez et al.||2013 > in SED data 
space. The top panel shows 30 model SEDs for both objects (gray 
solid lines), generated from a random sampling of the posterior 
parameter distributions shown in FigureOur photometric data 
is plotted on top (circles). The bottom panels show the residuals. 


between them. We can see that for each band, across all the 
posterior distributions, the differences between the models 
and the data never exceed ~ 0.1 mag. The discrepancies be¬ 
tween the model and data can be attributed to one or more 
of the following: inaccuracies in the intrinsic SEDs of OB 
stars in the Padova isochrones; inaccuracies in the shape of 
the reddening law; a calibration offset between the optical 
and NIR catalogues. 

Table [^compares the stellar parameters of the 21 known 
OB stars derived in this study with those from [Vargas] 
Alvarez et al.| ( |2Q13| ). Here Aq has been converted to Ay 


and the VPHAS+ g band magnitudes have been converted 
to V band using the Sloan to Johnson conversion from Lup- 
ton (2005|^for ease of comparison. We also note that our 
SED-derived log(Teff) v alues are compared to spectroscopic 


values where available (Vargas Alvarez et al. 


2013 


Rauw 


et al. 2007| . Otherwise effective temperatures are derived 


from spectral types according to the temperature scales of 


^ https://www.sdss3.org/dr8/algorithms/sdssUBVRITransform.php 



^[mags] 


Figure 8. The difference between stellar parameters fo und in this 


study and those found by Vargas Alvarez et al. ( 2013| . The solid 


line shows zero difference while the dashed line shows the median 
difference. 


Martins et al. (2005) and 

Zorec Briot 

(1991|). We restrict 

our comparison to the results in Vargas A 

Jvarez et al. (2013) 

based on the [Fitzpatrick & Massa (2007) extinction curves. 


Figure!^ plots the difference between the values derived 
in the two studies. It must be noted that star #584 has 
not been included in this analysis as extreme blending has 
substantially affected its photometry (see Figure[^and Table 

i- 


A significant difference is found between the trans¬ 
formed V band magnitudes in VPHAS+ and HST of ^ 


0.18 mag, such that VPHAS+ is brighter. Vargas Alvarez 


et al. (2013) compare their empirical B and V band mea¬ 


surements with those of Moffat et al. (1991) and Rauw et al. 


(2007) and find that those ground based measurements are 
also systematically bright er, by 0.18 and 0.15 mag, an d by 
0.22 and 0.12 respectively. Vargas Alvarez et al. (2013) sug¬ 
gest that the difference may be due to source blending fol¬ 
lowing on from the effects of atmospheric seeing. If this 
were the case we would expect to find objects in the most 
crowded/blended region of the cluster to be consistently 
more discrepant. As we do not see this effect we suspect 
a real calibration difference. Hur et al. (2015) have also un¬ 
covered a similar problem but find good agreement between 
their optical photometry and that of Rauw et al. (2007). If 
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Table 2. Table comparing the derived stellar parameters of objects with known spectral type from Vargas Alvarez et al. (20131 with 
the results in this study. The ID given corresponds to the numer ation given by ^^g^sA^y ^ez et al. ( 2013| . Most of the effective 
temperatures in the HST column were derived spectroscopically by Vargas Alvarez et al. ( |2013| and uncertainties were given. The rest 
have no provided uncertainty as they were estimated from their spectral types using the temperature scales from [Martins et al.| ( |2005| > 
and jZorec Briotj ( |1991| > . 


ID 

ST 

VPHAS+ 

Av 

HST 

VPHAS+ 

R-v 

HST 

log(TefF) 

VPHAS+ HST 

VPHAS+ 

HST 

V 

VPHAS+ 

HST 

137 

04 V 

7 47+0.04 
' •^'-0.04 

7.41 + 0.22 

4 r)4+0'05 
4.UO-0 05 

3.84 + 0.07 

4 ^4+0.05 
4.D9-Q 06 

4.633 + 0.004 

13.43lS:“ 

13.19 + 0.45 

15.496 ++0.056 

15.591 + 0.006 

178 

04 V-III((f)) 

‘:!4+0-04 

6.38 + 0.07 

4-o3lS;S6 

3.93 + 0.03 

4 fJ4+o.05 

4.629 + 0.002 

1 9 qq+0.64 

10 . 00 _Q g 2 

11.79 + 0.16 

14.385 ++0.055 

14.490 + 0.004 

395 

07.5V 

7C+0.07 

'°-0.08 

6.92 + 0.07 

4.08iS;S7 

3.77 + 0.03 

^•‘^^-0.07 

4.544 + 0.000 

1 2 91 +^•‘^9 

iz.yi_Q 7 g 

12.78 + 0.18 

15.688 ++0.056 

16.019 + 0.062 

505 

08.5V 

1 q+0.05 
'^•-^^-0.06 

6.36 + 0.14 

'J-O4-0 06 

3.71 + 0.06 

4 5Q+0-07 

4.531 + 0.006 

14 57+0.93 
i+.o^ -0 95 

13.29 + 0.30 

15.889 ++0.056 

16.094 + 0.005 

528 

08 V 

6 72+^-®® 

' ^-0.07 

6.97 + 0.14 

4i.uZ-o 06 

3.99 + 0.05 

4 Kfj+0.09 

4.DD-0 08 

4.544 + 0.005 

1 4 Q 4 + 1-10 
-*-"•"+-0.89 

12.55 + 0.30 

15.571 ++0.056 

15.841 + 0.005 

548 

04 V 

fi 44+0.05 

6.48 + 0.10 

^•^^-0.06 

3.76 + 0.04 

4 61 +0-0® 
^•"-‘--0.07 

4.633 + 0.002 

13 11 +0-81 
10 . 11-0 90 

13.19 + 0.23 

14.361 + +0.055 

14.522 + 0.002 

549 

B 1 V 


6.09 + 0.08 

4 14+0.09 
^•-‘-^-0.08 

4.01 + 0.04 

4 4tl+9-08 

4.398 + 0.000 

1 2 5 Q+"-®^ 
iz.oy_o 50 

11.68 + 0.19 

15.485 ++0.056 

15.562 + 0.005 

584 

08 V 

^.OU_o 04 

6.19 + 0.05 

2 01+0.04 

3.73 + 0.02 

4 . 66 I 0.05 

4.544 + 0.002 

15 24+0-^2 
iO.z ^_0 04 

12.94 + 0.12 

14.195 ++0.055 

15.442 + 0.004 

620 

B 1 V 

r yy+0.09 
' '-0.09 

5.77 + 0.08 

^•^^-0.08 

3.82 + 0.04 


4.398 + 0.000 

13.46lS;“ 

12.56 + 0.19 

16.007 ++0.057 

16.086 + 0.006 

640 

09.5V 

c OO+0.05 

'4-'J'^-0.07 

6.37 + 0.05 

4 07 + 0.07 
-0.06 

3.73 + 0.02 

^•‘^'-0.08 

4.505 + 0.002 


13.11 + 0.13 

16.065 ++0.057 

16.234 + 0.006 

704 

04 V 

D.Uc5_o 05 

6.27 + 0.29 

4 Q 4 +O.O 6 

'j.y^-0,06 

3.76 + 0.12 

^•'^-‘--0.07 

4.681 + 0.008 

12 91+"-'^'^ 
iz.yi_o 85 

14.26 + 0.63 

13.844 ++0.055 

14.059 + 0.002 

714 

03 V 

^•^-'--o.os 

6.08 + 0.11 

4 C 7 +O.O 6 
' -0.06 

3.73 + 0.05 

4 .OZ -0 07 

4.643 + 0.000 

14 29+"^®® 

i+-^y-0.82 

14.53 + 0.26 

14.642 + +0.055 

15.017 + 0.003 

722 

06 V 

y oi +0.05 
'•^■^-0.06 

7.23 + 0.04 

4 04+0.06 
'^•^^-0.05 

3.65 + 0.01 

4 61 +0-0® 
^•"-‘--0.08 

4.584 + 0.001 

1 2 7Q+0-84 
iZ. <y_0.94 

12.04 + 0.11 

14.944 + +0.055 

15.060 + 0.030 

738 

05.5V 

r 04 + 0 .05 

^•°4-0.05 

6.02 + 0.08 

Q OO+0.06 
'^•°°-0.06 

3.73 + 0.04 

4 61 +0-06 
^•'^-‘--0.07 

4.602 + 0.000 

1 3 91 +0-81 
ic5.yi_o 90 

13.39 + 0.19 

14.696 ++0.055 

14.896 + 0.003 

769 

09.5V 

D.OU_o 09 

6.63 + 0.06 

q C« + 0.07 

d.5D-o 06 

3.65 + 0.02 

4 54 + 0+0 

4.491 + 0.002 

14 12“^^'^" 

I+.IZ-O 89 

13.04 + 0.13 

16.351 ++0.057 

16.576 + 0.008 

804 

06 III 

7 11+0-11 

6.91 + 0.04 

4 n+0.09 
^•-‘--‘--0.09 

3.71 + 0.01 

^•°’0-0.08 

4.582 + 0.001 

12-14+J? 

11.78 + 0.10 

14.290 ++0.055 

14.433 + 0.003 

857 

04.5 V 

'^•‘^‘-’-0.08 

6.13 + 0.08 

4 1 7+0.09 
'-0.08 

3.63 + 0.03 

4 5fi+0-09 

+ -O°-0.08 

4.623 + 0.002 


12.65 + 0.18 

13.335 ++0.055 

13.869 + 0.003 

879 

09.5V 

c 77+0.06 
' '-0.07 

6.98 + 0.07 

q Qq + 0.06 
O.OZ-O 06 

3.70 + 0.03 

4 57 + 0.08 
' -0.08 

4.519 + 0.003 

14 37+1-05 

i+.o^ -0 98 

13.11 + 0.16 

16.510 ++0.058 

16.645 + 0.056 

913 

03-4V 

a 90 + 0.05 

u.zo-o 06 

6.42 + 0.11 

q 07 +O.O 6 
' -0.06 

3.66 + 0.04 

4 61 

+ •"-‘--0.07 

4.642 + 0.002 

1 3 9q+0.78 
io.zo_o g9 

13.45 + 0.24 

14.344 + +0.055 

14.531 + 0.002 

924 

08 V 

c 94+0.06 
^•^‘^-0.07 

6.40 + 0.07 

q CQ + 0.06 
'^•°°-0.06 

3.60 + 0.03 

4 57+0-08 
+ •"'-0.08 

4.544 + 0.000 

14.08lJ;“ 

13.16 + 0.16 

15.680 ++0.056 

15.960 + 0.005 

1039 

04-5 V 

44+0.05 

’^•^'^-0.05 

6.42 + 0.10 

q orj+0.06 

'^•°^-0.05 

3.47 + 0.04 

+ •"•^-0.07 

4.622 + 0.002 

13.261";™ 

12.98 + 0.22 

14.429 + +0.055 

14.523 + 0.030 


the scale of Rauw et al. (2007) is the right one, our photom¬ 
etry may be too bright by ~ 0.05 mag. 

The apparent systematic calibration difference between 
the two data sets is reflected in the derived values of Ay - In 
particular the median of the star-by-star differences in Ay 
shows that our e xtinctions are on average 0.14 mag less than 
those derived by Vargas Alvarez et al. (2013). The median 
Ay with the 16^^ and 84^^ percentiles in this study and in 


Vargas Alvarez et al. (2013) are; Av^(VPHAS+) = 6.34 


+0.44 


and A+(HST) — 6.41^Q;5g. With brighter optical magni¬ 
tudes there is also an offset in Ry such that our values are 
higher: the median star-by-star difference in Ry is 0.20, 
while sample medians are respectively Ra/(VPHAS+) = 


3.96;^;!^ and Ry(RST) = 3.73 


-0.08 
+ 0.11 • 


Despite the expectation of poor constraints on distance, 
the difference in the median values of fi happen to be very 
small: /i(VPHAS+) = and ^((HST) = 13.07:^J:g^ 

This is likely to be due to the O stars in Wd 2 being on the 
main sequence, matching our assumption. Similarly there is 
only a modest offset on average in the measures of effective 
temperature. 

The results of this comparison are encouraging. We have 
found good quantitative agreement, within the uncertain- 
ties, between our de rived parameters and those of j Varg^ 


Alvarez et al. 


(2013) drawing on HST optical photometry. 


Where there are differences, we understand their origin. This 
gives us confidence that both our method and the underlying 
VPHAS+ data are producing reliable results. 


5 RESULTS 

Here we apply the SED fitting methods discussed above to 
the full selection of OB candidates from our pilot ~2 sq.deg 
field. 


5.1 ‘Goodness-of-fit’ 

The posterior distributions obtained tell us the most prob¬ 
able parameters given the data, however they do not tell 
us anything about ‘goodness-of-fit’. As some objects in our 
selection may be contaminants or may just have bad pho¬ 
tometry, it is important to determine how well the data fit 
the model in order to obtain a ‘clean’ selection of OB stars. 
We have opted to use the value of given by the SED 
fits, at the median values in the marginalised posterior dis¬ 
tribution. We are aware that the posterior medians may not 
exactly trace the maximum likelihood, but they provide a 
representative sample. 

Eigure shows the distribution of the fits to all 
1050 objects in the wider selection above the distribution 
obtained for the known objects from Vargas Alvarez et al 


(2013). Since we are fitting 7 data points with 4 parame¬ 


ters we expect a /c = 3 x distribution peaking at 1 - the 
top panel of Eigure indicates this is what happens and, 
by implication, that the uncertainties on our data points 
are not significantly over- or under- estimated. In keeping 
with this, we have chosen to use the commonly adopted 5% 
significance level, at x^ = 7.82 as the limit beyond which 
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Figure 9. distributions for the known objects (bottom) and 
the wider selection (top). The x^ distribution for the wider selec¬ 
tion peaks at ~ 1 as a expected from a distribution with k = 3 
degrees of freedom. Using a 5% significance level we judge objects 
with x^ > 7.82 to be unsatisfactorily fit. The known objects with 
poor fits are subject to photometric blending in the cluster’s core. 


we judge the fits to the applied model to be unsatisfactory. 
This cut makes reasonable sense when applied to the x^ 


distribution for the known objects (in common with Vargas 
Alvarez et al.||2013 ), in that the 10 confirmed OB stars be¬ 


yond the chosen cut are mainly there because of the impact 
on the photometry of the blending in the crowded central 
parts of Wd 2 present in the VST data. For this reason we 
have still tabulated those objects not meeting our selection 
criteria but have not used them in any further analysis. We 


note that if both 2MASS Skrutskie et al. (2006) and Ascenso 


et al. (2007) photometry are available we keep which ever 
yields a better 


5.1.1 Further cross-matches with previously catalogued 
objects 

All of the objects in the initial selection were cross-matched 
to < 1" with the SIMBAD ( Wenger et al.|200'Q ) database to 
check for further examples of objects of known type. 


Tsujimoto et al. (2007) conducted a 17 x 17 arcmin high 


resolution X-ray imaging survey centred on Wd2 and the 
surrounding star forming region RCW 49. They identified 
17 new X-ray emitting O B candidates in this la r ger re gion, 
enclosing that studied by Vargas Alvarez et al. (2013). On 


using a 1” cross match radius we find 8 of these objects 
make it into our selection. Five of the missing objects are 
picked up by VPHAS+ but have ^ < 13 and hence were too 
bright to be selected. Conversely, the remaining 4 objects 
are detected by VPHAS+ but are too faint {g > 20) to be 
in our selection. It is likely that these objects are highly 
reddened. 

Across all other literature sources, accessed via SIM¬ 
BAD, fourteen further stars of confirmed type were found 
(see Table 1^. The breakdown of their classifications is as 
follows: six stars with a Wolf-Rayet (WR) component, three 
OV, two OIII, one OVb, one B5Vne, one carbon star and 


one star listed as MIIII. All six WR stars, the carbon star 
and one of the OV stars could not be fitted convincingly as 
reddened OB stars (i.e. x^ > 7.82), while the others were 
(x^ < 7.82). The OVb was confirmed a s an 0 3V + 05.5V 
binary system by Vargas Alvarez et al. (2013) but was not 
used in their SED fitting analysis - hence it did not feature in 
Section [4.2| On close inspection of the literature, it became 
clear that the SIMBAD MIIII attribution matching one of 
our selected objects is wrong, resulting from confusion over 
the sky position of the previously catalogued HAeBe can¬ 
didate, THA 35-11-41. THA 35-11-41 is indeed one of our 
selected objects but it is not at the position attributed to it 
by Carmona et al. (2010) where these authors observed an 
M giant spectrum. 

We also detect seven bright objects in the originally NIR 
selected open cluster DBS2003 45 ( Dutra et al.||200'3) cen- 


tred at 10hl9ml0.5s —58°02'22.6". The study by Zhu et al. 


(2009) identifies seven OB stars in this cluster estimated as 
ranging from spectral type BO to 07 from low resolution 
NIR spectroscopy. However, six out of seven of the positions 
given in Table 2 of Zhu et al. (2009) do not match with the 
VPHAS+ positions nor with any detections in the 2MASS 
point source catalogue. We therefore suspect that there is an 
error in the positions that they give whilst our objects are 
in common. We find these are among the most highly extin¬ 
guished objects in our selection with an average Ay = 8.37. 


5.1.2 Summary of results 

Figure shows the stages in the selection process: first, 
those stars without a match to good quality NIR photome¬ 
try have to be set aside (shown as grey crosses in the Fig¬ 
ure); next, those with ‘poor’ x^ values (the cyan-coloured 
squares); finally the good fits are divided in two groups based 
on their effective temperature. Those with a median poste¬ 
rior effective temperature exceeding 20000K, or equivalently 
log(Teff) ^ 4.3, are shown as red triangles while those that 
are assigned cooler fits are shown as blue squares. The hotter 
stars are our target group of spectral type B2 and earlier. 

Counter-intuitively perhaps, it can be seen in Figure 
that in the domain where g — r < 0.5, only 12 stars could 
be matched with good NIR photometry. This is because 
lowly reddened UV-excess objects detected in VPHAS+ 
are commonly too faint for detection in 2MASS due to 
their blue SEDs - for instance, some of these objects will 
be under-luminous hot compact objects. Unsurprisingly the 
cyan coloured squares representing objects with poor fits 
are frequently to be found above the Rayleigh-Jeans limit 
- only 2 objects with accepted fits fall into this part of the 
diagram. It is reassuring that there is some offset between 
the Rv = 3.1 B3V reddening vector, serving as lower bound 
to the selection region, and the spread of hotter objects: it 
suggests that few, if any, stars hotter than log(Teff) =4.3 
have been missed (given our other constraints, such as the 
magnitude limits). It is worth noting that the selection of 
objects that occupied the 0.1 mag wide band directly below 
the B3V reddening vector \n u — g provided just 1 star out 
of 374 with log(Teff) ^ 4.3 and x^ < 7.82. 

The main groupings emerging from the fitting process of 
all 1073 objects are shown in table[^ All of the objects along 
with their photometry and derived parameters are tabulated 
in Tables 111 and m 
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Table 3. Objects crossed matched with SIMBAD in the selection which have known spectral type. Derived parameters of highly evolved 
objects will be inaccurate due to the main-sequence assumption as shown by their large values. On further inspection of the literature 
the classification object ^895 is much different from that in SIMBAD (see Section ) 


ID 

RA 

DEC 

Identifier 

Spectral Type 

9 

log(Te//) 

Rv 

^0 

g 


282 

10 18 04.98 

-58 16 26.27 

WR 19 

WC5+09 

14.02 

4 07+0.05 

5 79+0-00 

' ^-0.09 

4 14+0.09 
^•-'■^-0.08 

9 60“*“®-^® 

y.DU-0 38 

39.60 

335 

10 18 53.39 

-58 07 52.94 

WR 19a 

WN 

15.45 


°-O^-0.09 

^•O^-0.06 

9 87“*“^-^® 

-0.89 

10.21 

437 

10 20 17.50 

-57 44 59.39 

C* 1665 

C* 

16.54 

^•OD_0.03 

lZ.CiU_Q Q4 

4.541^3 


436.47 

560 

10 22 05.75 

-57 53 46.03 

2MASS J10220574-5753460 

B5Vne 

15.71 

4 42+0-07 

5 fi4+0-09 

0-D4:_q gg 

3.80tHl 

11.98i°;g 

1.53 

644 

10 23 23.50 

-58 00 20.80 

SS 215 

02If*/WN5 

13.48 

4 41+0.06 

^•^-‘--0.05 

+0.09 
0-0 Wo.09 

4 27+0-10 
^•^'-0.09 

9 60“*“®-®^ 

y.DU-0 44 

16.06 

687 

10 23 58.01 

-57 45 48.93 

V* V712 Car 

03If*/WN6+03If*/WN6 

14.48 

^•^°-0.06 

^ .OU_o 09 

4 27+0-07 
Wo.07 

9 39“*“®-®^ 
y.oy-o 63 

10.47 

717 

10 24 01.20 

-57 45 31.03 

Cl* Westerlund 2 MSP 188 

03V+05.5V 

14.34 

4 tS4+0-lO 

a.oo_o 08 

6 7Q+0-09 

O- 'y_0.09 
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Figure 10. (u — g^g — r) diagram showing the stages of selection. 
Red triangles are the final selection used for further discussion. All 
of the objects clearly above the RJ reddening vector are returned 
as bad fits. 


5.1.3 Contaminants 


The (x^ > 7.82) fits have a range of causes. The most fre¬ 
quent are likely to be contact binaries or the products of 
poor photometry. 

Contact binaries may find their way into the selection 
because they are both quite common and rapidly variable. 
Figure 11 shows how around half of the x^ ^ 7.82 objects 


clearly separate in the (r — i, g — r) colour-colour diagram 
away from the OB stars towards redder g — r ai fixed r — i. 
This is plausibly the signature of contact binary (W UMa) 
interlopers. W UMa systems are doubly eclipsing binaries in 
which the brightness in any one band scarcely remains con¬ 


Table 4. Breakdown of the number of new OB candidates, previ¬ 
ously identified OB candidates and objects with known spectral 
type according to effective temperature and fit quality. 


All Objects: 1073 

log(Teff) ^ 4.3 

< 7.82 

x^ > 7.82 

Total 

527 

145 

New Candidate OBs 

489 

98 

Old Candidate OBs 

19 

28 

Known O - B2 stars 

19 

10 

Other 

0 1C star & 6 WR stars 

log(Teff) < 4.3 

X^ < 7.82 

x^ > 7.82 

Total 

321 

80 

New Candidate OBs 

320 

78 

Old Candidate OBs 

0 

2 

Known O - B2 stars 

1 

0 

All log(Teff) 

848 

225 


stant over time. These objects have typical orbital periods 


of 8 hours with two pronounced minima per cycle ( Rucin- 
ski]|1992 ). The u/g/r VST exposures are taken sequentially 


with about 15 minutes elapsing between u and g, and g and 
r. If the g band exposure of a W UMa system is taken at 
or near minimum light, its measured u — g colour is bluer 
than true, while g — r is redder, potentially pushing the star 
up into our OB selection. However these objects fail to pass 
as OB stars when the whole OnIR SED fit is performed, 
hence their poor x^ values. It has been estimated that there 
is around 1 W UMa system for every ~ 130 main sequence 
stars ( |Rucinski|1992| . So finding perhaps as many as ~ 100 
in our OB selection, given ~ 100000 stars across the 2 square 
degrees with ujgjr photometry, is reasonable. 

The second common origin for the poor fits is likely 
due to photometry affected by blending or incorrect cross¬ 
matching between bands. In the crowded core of Wd 2 this 
is an obvious difficulty (see Figures 4 and 10). 

The literature search already reported in section [5.1.2| 
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Ta ble 5: Sample table containin g the p ositions and photometry of all 1073 objects. The first five columns are the objects IDs given in this study, VPHAS ID, [Moffat et al.| ( |l99l] , MSP 
ID, Vargas Alvarez et al. (2013^, VA ID, Tsujimoto et al. (2007', TFT ID, and in SIMBAD, SIMBAD ID, where applicable. The full table can be found in the electronic version of this 
paper. 
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r — i 


Figure 11. Positions of objects with ^ 7.82 (blue dots) and 
> 7.82 (grey crosses) in the {r — i, g — r) plane. The solid 
black line is the reddening vector of an 09V with Ry = 3.8. 
The dashed line is the unreddened main sequence. We find that a 
large number of objects with ‘poor’ fits fall away from the OB star 
reddening vector. These objects show colours that are consistent 
with eclipsing W UMa contact binaries. 


revealed that high may be linked to extreme objects 
like WR stars (6 examples) and carbon stars (1 only). An¬ 
other rare contaminant may be white dwarf/M dwarf bina¬ 
ries that can present blue u — g, alongside red r — i. The 
blue white-dwarf light begins to be overwhelmed by the 
red dwarf’s light with increasing wavelength, shifting the 
combined colours below and to the right of the OB red¬ 
dening track in the [g — r^r — i) diagram (fig EH- Such 
objects are known to co-locate with reddened OB stars in 
the {u — g^g — r) diagram or they may fall beyond the RJ 
reddening vector ( Smolcic et al.||2004 ). 


distance modulus range and show a fairly wide spread in 
extinction law with 3.5 ^ Ry ^4.5. 

Echoing the initial mass function (IMF), the distribu¬ 
tion in median log(Teff) values is heavily skewed towards the 
lower end. The turn over in the log(Teff) distribution at just 
below log(Teff) = 4.3 further supports the conclusion that 
our initial selection of VPHAS+ sources in the (u — g,g — r) 
diagram is essentially complete in the desired O to B2 effec¬ 
tive temperature range (given our magnitude limits). The 
coolest object in the candidate list is ~ 16000 K. 

Stars with median estimated effective temperatures in 
excess of 30000K (log(Teff) ^ 4.477) are regarded as can¬ 
didate O stars. Of the new discoveries 74 meet this crite¬ 
rion. We can further subdivide this group to distinguish the 
highly probable O stars: 28 objects have a 16^^ percentile 
log(Teff) exceeding 4.477. Seven of these may be sdO stars 
(see section 5.3). 

Predictably, many of the hottest candidates are in and 
around Wd 2: this young massive cluster does indeed stand 
out in this part of the Galactic Plane. Moreover the top 
left panel in \T^ suggests a relative lack of cooler OB stars 
within the 8 arcmin box centred on the cluster. This could 
be taken to imply that the stellar mass function of Wd 2 and 
environs is top heavy. At the same time, there are biases that 
can favour the detection of more massive stars at the likely 
distance of the cluster (// ~ 13 — 14) - namely, the effects 
of crowding (less massive fainter stars are more likely to be 
lost in blends) and of magnitude limited selection. However 
this is unlikely to be all of the explanation given that there 
are plenty of examples of Aq ~ 5.5 — 7 cool candidates with 
a similar estimated distance modulus. 

Figure shows the upper and lower uncertainties on 
each parameter as a function of ^-band magnitude for all 
< 7.82 objects. The uncertainty on log(Teff) and Aq in¬ 
creases for fainter objects, tracking the increase with rising 
magnitude of the photometric errors. Conversely the uncer¬ 
tainty on Ry shows a slight increase with decreasing magni¬ 
tude at the bright end. Ry is more difficult to determine for 
bright objects as they tend to be less obscured. Neverthe¬ 
less it is evident that both Ry and Aq are consistently well 
determined across the entire magnitude range. Our OnIR 
SED fits deliver Aq to within < 0.09 mag up to 18^^ mag¬ 
nitude, rising up to < 0.25 mag at 20^^ magnitude. We find 
the median uncertainty on Ry to be 0.081. 


5.2 Parameters of the candidate OB stars 


Figure shows the distribution of stellar parameters across 
the entire selection for the objects fitting successfully to a 
reddened OB-star SED (x^ ^ 7.82 and log(Teff) ^ 4.3). 
Coloured in red are the results for all objects within an 8 
arcmin box centred on Wd 2 (drawn in Figure p~7| . It can be 
seen that those objects in or near the cluster are reported 
to have similar extinction in the range 5.5 ^ Aq ^ 7 (top 
right panel in Figure 12). Otherwise, the reddenings range 


more broadly across the full 2 square degrees from Aq ^ 3 
up to Ao ^ 8. Other features of this particular sight-line are 
that larger than standard Ry is favoured - a roughly normal 
distribution in Ry about a mean value of Ry = 3.84zb0.25 is 
obtained - and that most of the selected stars are attributed 
distances of between ~ 2 kpc (/x ~ 11) and ~ 6 kpc (/x ~ 14). 
The objects in/near Wd2 tend toward the higher end of the 


5.3 Inferences from the best-fit parameters and 
other aspects of the photometry 

A richer understanding of the candidate objects can be ob¬ 
tained from a combination of more scrutiny of the fit param¬ 
eters obtained and from a fuller utilisation of the VPHAS+ 
photometry at our disposal. So far the focus has been on 
the information to be extracted from the individual OnIR 
SEDs - treating all candidates as if they are well described 
as reddened, single, main-sequence OB stars. We can learn 
more through consideration of the ensemble of objects, and 
if use is made of the narrowband Ho band to separate out 
emission fine stars. 

First we acknowledge and relax the main sequence as¬ 
sumption applied so far. The first two panels of Figure 
show scatter plots of the best-fit median distance modulus, 
/X, vs. log(Teff) and vs. Aq for the candidate OB stars. Differ- 
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Candidates 

O OBs # Sub-luminous # Blue supergiants 


Known 

A OV ▼ OIII □ OV+OV ★ WR # BV • HAeBe 



4.3 4.4 4.5 4.6 4.7 4 6 8 10 0.0 0.8 1.6 2.4 3.2 


^og{Teff) ^0 g-r 

Figure 14. 2-D distribution of the best fit parameters for the final selection of OB candidates (x^ ^ 7.82 and log(Teff) ^ 4.3) and objects 
in the selection with known spectral type (the carbon star lays outside the range of 2 of these diagrams and was therefore not included). 
Objects shown in red and yellow are thought to be candidate sub-luminous OB stars and candidate blue supergiants respectively. The 
areas shaded in grey are where we cannot detect OB stars given the survey limits. 



Rv fji 

Figure 12. Distribution of the best fit parameters for the selec¬ 
tion of objects with < 7.82. The red bars are objects within 
and 8 arcmin box of Wd 2 while the grey bars are the wider se¬ 
lection. We find that the objects spatially associated with Wd 2 
show a tight distribution in Aq and provide an over density of 
objects in the 5.5 ^ Aq ^ 7 range and also show a wider spread 
in Ry. 
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Figure 13. Uncertainty on each parameter as a function of g 
band magnitude. Uncertainties are derived from the 16^^ and 
84 th percentiles of the posterior distributions. 


ent symbols are over-plotted to pick out the already known 
objects listed in SIM BAD as well as the O stars of | Vargas] 


Alvarez et al. 


(2013). The areas shaded in grey are where we 


cannot detect OB stars given the survey limits. The objects 
plotted as red circles have relatively low extinction but, if 


we take the returned distance moduli at face value, they 
would have to be construed as very distant (> 10 kpc) when 
compared to the known OB stars. It is more plausible that 
these are intrinsically sub-luminous objects rather than dis¬ 
tant OB stars located in remarkably clear reddening holes. 
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Their scattered spatial distribution across the whole field 
shown in Figure pT| supports this argument. 

The converse argument can be applied to those objects 
plotted as yellow circles: they are found to have more than 6 
magnitudes of extinction but are seemingly very close (less 
than ~700pc away, /x < 9). We suspect that these objects 
are intrinsically much higher-luminosity, evolved B stars. 
The proximity of these stars in the figures to the (poorly- 
fit) known WR stars, including the highly-luminous WR20a, 
lends credibility to this interpretation. 

Referring back to the photometry in the form of a 
{g^g — r) colour magnitude diagram (CMD), these inter¬ 
pretations are seen to make sense - the sub-luminous and 
over-luminous objects form tracks separated from the main- 
sequence - see the third panel of Figure Table lists 
these extreme objects. There will be further discussion of 
them in Section 6. The main concentration of objects ap¬ 
pears in the 11.5 < g < 14 mag range which equates to 
distances ranging from 2 — 6kpc. This encloses the derived 
distance range of the Carina arm traced in CO by |Grabelsky| 


et al. 


(1988), near its tangent. 


We can also use the VPHAS+ Ho measurements to un¬ 
cover any emission line stars in our selection. The presence 
of emission lines implies the presence of ionized circumstellar 
gas which, among massive OB stars, most commonly indi¬ 
cates classical Be stars with circumstellar disks. Although 
the OnIR SEDs of classical Be stars are not greatly differ¬ 
ent from normal B stars of similar effective temperature, 
the derived interstellar extinctions from SED fits that do 
not take into account the circumstellar continuum emis¬ 
sion will nevertheless be overestimated. We have used the 
(r — i, r — Ho) diagram to select all objects that lie more 
than 0.1 mag in r —Ho above the 09V reddening vector (this 
equates to ~ 10A in emission line equivalent width). Eigure 


15 shows this selection. Using the relation between EWiYia) 


and added colour excess E{B — V) due to the presence of 
a circumstellar disk in classical Be stars from IDachs et ah] 
(1988), we can estimate that the derived reddenings (Ao) for 
our Ho-excess stars will have been inflated by between ~0.1 
and ~0.3 magnitudes. There are 17 of these objects in the 
^ 7.82 and log(Teff) ^ 4.3 group and a further 63 with 
X^ > 7.82 and/or log(Teff) < 4.3. Objects with Ho excess 
are marked in Tabled 


5.4 Reddening 

After removing the obvious sub/over-luminous objects and 
the emission line stars from the selection we are left with a 
cleaner selection of 458 ~ non-emission OB candidates and 
19 known OB stars available for further examination of their 
reddening properties. 

Given our tight grasp on Ao and Rv, it is of interest to 
consider their interdependence. Rv is plotted as a function 
of Ao in Eigure The left panel of this Eigure includes 
those objects within an 8 arcmin box around Wd 2 and 
the right hand panel excludes them. The areas shaded in 
grey are where we cannot detect OB stars given the survey 
limits. In both cases we can see a moderate positive cor¬ 
relation in Rv as a function of Ao (correlation coefficient 
r = 0.47 and r = 0.45 respectively). On comparing the two 
panels, it is evident that the members of Wd 2 drive up the 
Rv trend more sharply when they are included. The shaded 


0.9 


0.6 


0.0 


Figure 15. 13 of the x^ ^ 7.82 and log(Teff ) ^ 4.3 objects show 
Ho; excess. As emission is usually associated with circumstellar 
dust; the derived extinction may be incorrect. The solid line is 
the reddening vector of an 09V raised by 0.1 in r — Ha. 
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background shows that the trends seen are independent of 
the boundaries set by the survey selection limits. Given that 
it was demonstrated in Section [4^ that the fitting method 
generates negligible covariance between Ao and Rv, we can 
say with confidence that the correlation apparent now is re¬ 
lated to the physical nature of the volume of space under 
study. 

It is commonly accepted that increasing Rv is linked 
to increasing typical dust grain size, and that values of 3.5 
and more are associated with denser molecular cloud envi¬ 
ronments (see e.g. Draine|2003 ). The ~2 square degrees un¬ 
der examination here sample sight-lines lying just inside the 
Garina Arm tangent direction. Our pencil beam is evidently 
one that would initially pass through the atomic diffuse in¬ 
terstellar medium and then enter the dense clouds of the 
Garina Arm, wherein Wd 2 is located. In this situation it 
makes sense that as the dust column grows it becomes ever 
more dominated by the dense/molecular ISM component - 
i.e. Rv tends to rise. However the rise is not dramatic, and 
the data points show significant dispersion, which may imply 
that the variation in the dust properties within the sampled 
volume is not especially coherent. The effect of the bright 
limit of the survey is to remove sensitivity to Ao much below 
2-3, or to distances less than ~3 kpc (see below). Gurrent 
maps of Galactic spiral arm structure place this distance 
already within the Garina Arm ( Russeil|2003| Vallee|2014 ). 

The clear message of Eigure |l6| is that the typical, 
if necessarily idealised, reddening law for this sight-line is 
Rv ~ 3.8, which rises much less sharply with decreasing 
wavelength than the Galactic average of Rv = 3.1 (see Eig¬ 
ure 13 in Eitzpatrick Massa|2007 ). 
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Table 7. Table containing the derived stellar parameters of the sub-luminous and blue supergiant candidates in the ^ 7.82 and 
log(Teff) ^ 4.3 group. 


ID 


RA 


DEC 

g 

log(Teff) 

^0 

Rv 



Sub-Luminous 

19 

10 

13 

40.87 

-57 

43 

15.38 

18.13 

4 65+°-°^ 

^.DO_o 05 

4 32“^®'®® 

4 95“^®'^® 

^•^^-0.15 

18.82+;“ 

7.09 

28 

10 

13 

54.74 

-58 

14 

33.73 

18.91 

4 oi+0.05 

4 30+0-10 

4.OD-0 13 

q OO + 0.17 
'^•°^-0.16 

15.441?;^? 

3.53 

82 

10 

15 

21.46 

-57 

53 

30.64 

19.43 

4 07+0.08 
-0.06 

5 94+0.13 
^•^^-0.13 

4:.UO_o 10 

14.82+;^^ 

1.11 

86 

10 

15 

22.76 

-57 

48 

21.73 

18.67 

-0.08 


4 17+0-23 
^•1' -0.21 

1 8 qq+0.98 

10.00_Q gg 

2.10 

89 

10 
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With Wd2 


Without Wd2 





Figure 16. Aq vs. Ry plot for the final selection with and without the objects within the 8 arcmin box surrounding Wd 2. There modest 
increase in Ry as a function of Aq in either case with correlation coefficient r = 0.47 and r = 0.45 respectively. The areas shaded in grey 
are where we cannot detect OB stars given the survey limits. 


6 DISCUSSION 

6.1 The number and spatial distribution of the 
OB candidates 

Figure shows the location of each new candidate in the 2 
square degrees for which the SED fit returned ^ 7.82 and 
log(Teff) ^ 4.3, over-plotted on the VPHAS+ Ho mosaic. 
Each star is colour-coded according to its derived extinction, 
Aq. The 527 objects are scattered across the field, with lower 
reddenings (Aq < 5) dominating in the southern half. Apart 
from in Westerlund 2 itself, the distribution is sparser and 
more highly-reddened in the north. Towards the NW and 
the tangent direction, roughly at RA lOhllm, Dec -56 14 
(J2000) ( |Dame IpOOTt , the most reddened objects (Aq > 8) 
are found. 

489 of the objects shown have not been identified previ¬ 
ously as confirmed or candidate OB stars. Previous works by 


Reed (2003), and by Kaltcheva & Golev (2012) have noted 
a further 26 stars earlier than B3 within this region - all of 
which are brighter than V =11, and therefore not in our 
sample. Also for reasons of brightness, our sample does not 
include any stars obviously associated with IC 2581. Turner 
(1978) studied the cluster - home to a number of early B 
stars - establishing a distance of 2.87 kpc, and a typical red¬ 
dening corresponding to Aq ~ 1.5. Eor the present selection, 
this cluster is too close and too lowly-reddened: a B3 main 
sequence star with Aq = 1.5 needs to be at a distance of 3.8 
kpc to achieve ^ = 13. The one star that has been uncov¬ 
ered close to IC 2581 is a candidate sub-luminous object, 
likely to be at a much shorter distance and unconnected to 
IC 2581. It can be seen in Eigurep^that Aq = 1.9 for the 
least reddened candidate OB star in the sample. 


It has been argued before by e.g. Grabelsky et al. (1988) 


and Dame (2007) that the Carina Arm tangent region traced 


the Solar Circle where declining amounts of molecular gas 
are detected. Taking note of these considerations, it seems 
likely that the high values of Ry, trending from 3.6 to 3.9, 
revealed by our SED fits (Eigure |16[ ) are largely a product 
of the dominant and increasing contribution to the total ex¬ 
tinction from the dust column of the Carina Arm. Similarly 
Povich et al. (2011) found it necessary to adopt Ry = 4 
for embedded Carina Arm objects at I ~ 287°. In contrast. 
Turner (1978) determined Ry to be 3.11zb0.18 across this 


region, based on bright OB stars with extinctions below 
Ay of 2 - clearly the foreground to our sample. Indeed it 
seems likely that much of the extinction of the OB popu¬ 
lation spanning Aq ~ 2 to Aq ~ 9 accumulates within the 
Carina Arm. The appearance of Eigure indicates few de¬ 
tected main sequence OB stars beyond a distance modulus 
of 14 ( 6kpc). 


6.2 Westerlund 2 

Eigures p^ and p^ tell us that a single value of Ry cannot be 
used to describe the extinction law of sight lines towards all 
objects in Wd 2. Instead we find that Ry ranges from ap¬ 
proximately 3.5 to 4.5 within the cluster. Similar spreads in 
Ry within star clusters has previously been found by |Eitz-| 
Patrick & Massa (2007) and highlights the importance of 
deriving Ry on a star-by-star basis. Hur et al. (2015) de¬ 


in CO spans the distance range from 3 to 5 kpc. At larger 
distances the conical volume captured here reaches beyond 


scribe a hybrid extinction model with Ry = 3.33 zb 0.03 to 
Ao ~ 3 (based on three stars), while Ry = 4.14 zb 0.08 is 
required for stars in Wd 2. Eigure p^ cautions against this 
clear-cut interpretation, even while our results are numer¬ 
ically consistent with theirs. Reality is more fractal and it 
is best not to place too much weight on small numbers of 
stars. 

In Eigure p^ we noticed a tight distribution in Aq for 
the objects close to Wd 2 as projected on the sky. While 
there are no new OB star candidates in the central region of 
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Table 8. The reddening parameters and angular separation from 
the centre of Wd 2 ( RA 10 24 18.5 DEC -57 45 32.3 (J2000)) 
of new O star candidates with similar reddening to the cluster, 
outside the 8 arcmin box shown in figure |17| All objects have 
log(Teff) > 4.477 and 5.8 > Aq > 7.2. See Tablesandfor the 
full set of data. 


ID 

9 

Ao 

Separation 

(arcmin) 

44 

16.62 

6-881°:°^ 

80.21 

121 

14.28 

6-2llao8 

70.59 

191 

17.43 

D.UO_o 09 

64.98 

144 

16.94 

« 1 O+0.07 
'^•-^'^-0.09 

68.46 

161 

19.28 

71-1 -|-0.07 

62.73 

346 

16.60 

' °-0.09 

46.51 

413 

17.06 

Qq+0.05 

b.8C)_o 06 

36.42 

576 

15.40 

6 42+0-05 

O-^^_0.07 

23.13 

916 

15.01 


19.67 

646 

15.58 

O-O^-0.08 

12.60 

796 

16.37 

7 nq+0.05 

‘ •’O'J-0.07 

12.50 

662 

15.33 

D.<:)U_o 00 

12.06 

846 

16.82 

O.Ob»_o 05 

6.03 

661 

17.21 

Ql +0.04 
^•^-*--0.05 

5.05 


Wd 2 (within the 8 arcmin box shown in fig 17), there are 
a handful of probable O stars scattered across the field that 
share its extinction. All objects that have extinctions con¬ 
sistent to within 1 cr of the mean of known objects in Wd 2 
(5.8 > Ao > 7.2) and have log(Teff) > 4.477 are identified in 
Table It is possible that these have been ejected from Wd 
2 by dynamical interactions or after supernova explosions in 
binary systems ( Allen Poveda|1971 Gies Bolton|19^ ). 
Given a derived distance of ~ 5 kpc to Wd 2, an object that 
is separated from the cluster by ~ 20 arcmin on the sky 
would have to have travelled a minimum distance of ~ 30 
pc in 1-2 Myr. This would equate to a minimum (plane-of- 
sky) velocity of ~ 25 km/s. Given that massive stars can 
attain runaway velocities of up to ~ 200 km/s through dy¬ 
namical encounters between binary systems ( [Gvaramadze 


et al. 2010), it is not unreasonable to consider that these 


objects have been ejected recently from Wd 2. Alternatively 
these stars may have formed in situ within the wider star 
forming region on a similar time scale to the cluster. Low 
resolution spectroscopy of object ^916 and ^646 suggests 
their spectral types may be early as 03-04 (Mohr-Smith et 


al in prep). Their positions are marked on Figure 17 


As a final point of interest, we note that the more red¬ 
dened cluster DBS2003 45, picked out in figure pT| as well, 
also appears to be surrounded by a scatter of OB stars that 
are reddened similarly to the cluster. 


6.3 Candidate blue supergiants and sub-luminous 
stars 


The results from Section [5i3| suggest the presence of 5 high 
luminosity B stars scattered across the field. If they are 
early-B supergiants, their absolute visual magnitudes would 
be in the region of ~ —6.5 (Growther et al. 2006). On 


correcting the previous main-sequence assumption, we find 
their derived distance moduli, /i, rise from ~ 9 to ~ 13.5, 
placing them amongst the general OB population that we 
pick out. Meylan & Maeder (1983) estimate a surface density 
of around 10 - 20 blue supergiants (BSGs) per kpc^ in the 
Galactic Plane. Assuming that our selection spans distances 
from 2-6 kpc we are sampling a projected disk surface area 
of a little over 1 kpc^; so finding 5 candidates undershoots 
the surface density prediction but not to the extent that it 
can be claimed to be inconsistent with it. Given that these 
candidates are affected by saturation in the i band {i < 12), 
there may one or two BSGs that have fallen into the ‘poor- 
fit’ group due to saturation in one or more bands. 

We also find evidence for the presence of a population of 
subdwarf stars (see table and figure p^. Of these 9 may be 
sdO stars, leaving 23 in the sdB category. The absolute mag¬ 


nitudes of the latter range from My = 3 — 6 (Stark & Wade 


2003). Since these objects are ~ 6 mag fainter than their 


main-sequence counterparts, their distance moduli are likely 
to be ~ 10 as opposed to the estimated values of ~ 16. This 
behaviour and the spatial scattering of the subdwarf candi¬ 
dates suggests that we are looking at a group of moderately 
reddened Aq ~ 4 stars in the foreground of the main OB 
population. We are biased to select more highly reddened 
subdwarf stars due to the 2MASS faint limit as discussed in 
Section |5.1.2| If this limit was not in place we would expect 
to find more lowly reddened sdB stars in the selection. 

Although the SED-fitting we have performed has no 
sensitivity to surface gravities and limited sensitivity to stel¬ 
lar effective temperature, the fact that the Garina Arm re¬ 
gion studied falls near the tangent has allowed us to pick 
out the extreme objects purely from their outlying distance 
moduli - relative to the near MS stars concentrated in the 
range 11 < /x < 14. While this approach works here, it is ev¬ 
ident that in other sight-lines, where the population of OB 
stars may be spread more uniformly across a larger distance 
range, the luminosity extremes would not stand out in the 
same way. 


7 CONCLUSIONS 

In summary, we have demonstrated a method for selecting 
and parametrizing the reddening and basic stellar properties 
of OB stars uniformly across large areas of the Southern 
sky using VPHAS+, and NIR survey data. The selection 
presented here has resulted in reddening parameters for 848 
O and B stars (see table[^. Of these, 489 are well-fit new OB 
stars hotter than 20000 K, including 74 probable O stars. 
This has been achieved by reaching down to ^ = 20 mag 
and approaches a factor of 10 increase relative to the small 
number of known and candidate O to B2 stars in the region. 

By bringing together VPHAS+ u,g,r,i photometry 
with NIR 2MASS photometry, we are able to determine 
both the value of the extinction, Ao, and test and select the 
most appropriate reddening law, as parametrised by Rv, to 
a high degree of accuracy: both are typically measured to 
better than 0.1 (magnitudes in the case of Ao). Pleasingly 
there are signs that the still preliminary nature of the photo¬ 
metric calibration of the VPHAS+ survey data blends well 
with the now well-established 2MASS calibration. 

We set out expecting to only gain a crude impression of 
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stellar effective temperatures (and hence distance moduli), 
and so it has turned out. But we have found a satisfying con¬ 
sistency with earlier results in our benchmark region around 
the much-studied cluster Westerlund 2, confirming that our 
methods are sound and able to e.g distinguish early O stars 
from late-0 and early-B stars. This represents an efficient 
start to selection that needs to be followed up by spectro¬ 
scopic confirmation and measurement of stellar parameters. 
With precise spectroscopic parameters in hand, the photom¬ 
etry can be re-used for direct and even more precise mea¬ 
surement of reddening laws. 

We have also seen how the high resolution and wide 
field of view of OmegaCam can bring a wider context to 
the study of open clusters and OB associations, through an 
ability to identify potentially-related stars that have either 
been ejected from clusters or simply have formed - perhaps 
as part of a wider star-format ion event - in relative isolation 
in the surrounding held. This study has also uncovered 5 
BSG candidates as well as 32 reddened candidate subdwarfs 
of which 9 may be sdO stars. 

In the future, we aim to roll out this method to support 
the complete characterisation of the massive-star population 
and the patterns of extinction they can reveal across the 
entire Southern Galactic mid-plane to distances of ~ 5kpc 
or more. Garmany et al. (1982) were able to claim a volume- 
limited census to ~ 2.5kpc 3 decades ago - now it should 
be possible to expand the effective volume by a factor of 4 
or so, with the difference this time that Gaia parallaxes as 
they appear will bestow a conhdence as to what the volume 
limits actually are. 


ACKNOWLEDGEMENTS 


This paper is based on data products from observations 
made with ESO Telescopes at the La Silla Paranal Observa¬ 
tory under programme ID 177.D-3023, as part of the VST 
Photometric H Survey of the Southern Galactic Plane and 
Bulge (VPHAS+, www.vphas.eu). Use is also made of data 
products from the Two Micron All Sky Survey, which is 
a joint project of the University of Massachusetts and the 
Infrared Processing and Analysis Genter/Galifornia Insti¬ 
tute of Technology, funded by the National Aeronautics and 
Space Administration and the National Science Foundation. 

This research made use of Astropy, a community- 


developed core Python package for Astronomy (Astropy 
Gollaboration et al.||20I3 ) and TopGat ( Taylor||2Q05 ). 

The authors wish to thank the referee for helpful 
comments that have improved this paper. MM-S acknowl¬ 
edges a studentship funded by the Science and Technol¬ 
ogy Facilities Gouncil (STFG) of the United Kingdom (ref. 
ST/K502029/I). JED and GB acknowledge the support of 
a research grant funded by the STFG (ref. ST/J00I333/I). 
NJW is in receipt of a fellowship funded by the Royal As¬ 
tronomical Society of the United Kingdom. 


REFERENCES 

Allen C., Poveda A., 1971, Astrophysics and Space Science, 13, 
350 

Ascenso J., Alves J., Beletsky Y., Lago M. T. V. T., 2007, A&A, 
466, 137 


Astropy Collaboration et al., 2013, A&A, 558, A33 
Bertelli G., Bressan A., Chiosi C., Fagotto F., Nasi E., 1994, 
A&AS, 106, 275 

Bestenlehner J. M. et al., 2011, A&A, 530, LI4 
Bressan A., Marigo P., Girardi L., Salasnich B., Dal Cero C., 
Rubele S., Nanni A., 2012, MNRAS, 427, 127 
Cardelli J. A., Clayton G. C., Mathis J. S., 1989, ApJ, 345, 245 
Carmona A., van den Ancker M. E., Audard M., Henning T., 
Setiawan J., Rodmann J., 2010, A&A, 517, A67 
Carraro G., Turner D., Majaess D., Baume G., 2012, ArXiv e- 
prints 

Crowther P. A., Lennon D. J., Walborn N. R., 2006, A&A, 446, 
279 

Dachs J., Kiehling R., Engels D., 1988, A&A, 194, 167 
Dame T. M., 2007, ApJL, 665, LI63 

de Bruijne J. H. J., 2012, Astrophysics and Space Science, 341, 
31 

de Wit W. J., Testi L., Palla F., Zinnecker H., 2005, A&A, 437, 
247 

Draine B. T., 2003, ARA&A, 41, 241 
Drew J. E. et al., 2014, MNRAS, 440, 2036 
Dutra C. M., Bica E., Soares J., Barbuy B., 2003, A&A, 400, 
533 

Fitzpatrick E. L., Massa D., 2007, ApJ, 663, 320 
Foreman-Mackey D., Hogg D. W., Lang D., Goodman J., 2013, 
PASP, 125, 306 

Garmany C. D., Conti P. S., Chiosi C., 1982, ApJ, 263, 777 
Gies D. R., Bolton C. T., 1986, ApJS, 61, 419 
Grabelsky D. A., Cohen R. S., Bronfman L., Thaddeus P., 1988, 
ApJ, 331, I8I 

Gvaramadze V. V., Gualandris A., Portegies Zwart S., 2010, 
in lAU Symposium, Vol. 266, lAU Symposium, de Grijs R., 
Lepine J. R. D., eds., pp. 4I3-4I6 
Hur H., Park B.-G., Sung H., Bessell M. S., Lim B., Chun M.-Y., 
Sohn S. T., 2015, MNRAS, 446, 3797 
Johnson H. L., Morgan W. W., 1953a, ApJ, II7, 313 
Johnson H. L., Morgan W. W., 1953b, ApJ, II7, 313 
Kaltcheva N. T., Golev V. K., 2012, PASP, 124, 128 
Langer N., 2012, ARA&A, 50, 107 

Maiz-Apellaniz J., Walborn N. R., Galue H. A., Wei L. H., 2004, 
ApJS, I5I, 103 

Martins F., Schaerer D., Hillier D. J., 2005, A&A, 436, 1049 
Meylan G., Maeder A., 1983, A&A, 124, 84 
Moffat A. F. J., Shara M. M., Potter M., 1991, AJ, 102, 642 
Parker R. J., Goodwin S. R, 2007, MNRAS, 380, 1271 
Portegies Zwart S. F., McMillan S. L. W., Gieles M., 2010, 
ARA&A, 48, 431 

Povich M. S. et al., 2011, ApJS, 194, 6 
Rahman M., Murray N., 2010, ApJ, 719, 1104 
Rauw G., Manfroid J., Gosset E., Naze Y., Sana H., De Becker 
M., Foellmi C., Moffat A. F. J., 2007, A&A, 463, 981 
Rauw G., Sana H., Naze Y., 2011, A&A, 535, A40 
Reed B. C., 2000, AJ, 120, 314 
Reed B. C., 2003, AJ, 125, 2531 
Rucinski S. M., 1992, AJ, 104, 1968 
Russeil D., 2003, A&A, 397, 133 
Skrutskie M. F. et al., 2006, AJ, I3I, 1163 
Smartt S. J., 2009, ARA&A, 47, 63 
Smolcic V. et al., 2004, ApJL, 615, LI4I 
Stark A. A., Lee Y., 2005, ApJL, 619, LI59 
Stark M. A., Wade R. A., 2003, AJ, 126, 1455 
Taylor M. B., 2005, in Astronomical Society of the Pacific Con¬ 
ference Series, Vol. 347, Astronomical Data Analysis Software 
and Systems XIV, Shopbell P., Britton M., Ebert R., eds., p. 29 
Tsujimoto M. et al., 2007, ApJ, 665, 719 
Turner D. G., 1978, AJ, 83, 1081 
Vallee J. R, 2008, AJ, 135, 1301 
Vallee J. R, 2014, AJ, 148, 5 













New OB stars in the Carina Arm 21 


Vargas Alvarez C. A., Kobulnicky H. A., Bradley D. R., Kan- 
nappan S. J., Norris M. A., Cool R. J., Miller B. P., 2013, AJ, 
145, 125 

Wenger M. et al., 2000, A&AS, 143, 9 

Zhu Q., Davies B., Figer D. F., Trombley C., 2009, ApJ, 702, 
929 

Zinnecker H., Yorke H. W., 2007, ARA&A, 45, 481 
Zorec J., Briot D., 1991, A&A, 245, 150 



